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Abstract 

We find that circular kinks form on the surface of granular material when the 
axis of rotation is tilted more than the angle of internal friction of the material. 
Radius of the kinks is measured as a function of the spinning speed and the 
tilting angle. Stability consideration of the surface results in an explanation 
that the kink is a boundary between the inner unstable and outer stable 
regions. A simple cellular automata model also displays kinks at the stability 
boundary. 
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Granular materials behave differently from any other familiar forms of matter. They 
possess both solid- and fluid-like nature and exhibit unusual dynamic behaviors, such as 
segregation, surface waves, heap formation and convection The surface of granular 

material in a spinning bucket is an example of such interesting phenomena. Vavrek and 
Baxter |^ showed that the surface shape of sand in a vertical spinning bucket can largely be 
explained using Coulomb's criterion. Medina et al. [0 investigated hysteresis of the surface 
shape with semi two-dimensional rotating bins and showed the existence of multiple steady- 
state solutions. Yeung studied the system with an initially conical surface. By using a 
model of granular surface flow |]9|JT0[|, he found that the behaviors of the model agreed well 



with the experiments. 

When we tilt the rotational axis of a bucket from the vertical direction, circular kinks 
develop on the granular surface if the tilting angle is greater than the angle of internal 
friction of the material. A glass beaker of 10.5 cm diameter and 1 liter capacity is used 
as a spinning bucket. The beaker is mounted on a dc motor to rotate around its axis of 
symmetry. The motor is fixed to a stand such that the tilting angle a can be varied. The 
rotation rate is varied from to 300 rpm by controlling the voltage. The angular velocity uj 
is measured by using a photogate timer and is constant within an error smaller than 1 rpm. 
We use natural sand, as used in general construction, as a prototypical granular material. 
In order to insure the monodispersity of sands, two sieves with 0.35 and 0.25 mm meshes are 
used and the sizes in between are selected. The angle of internal friction (angle of repose) 
and apparent density are found to be ^/ = 34 ± 1° and p = 1.52 g/cm^, respectively. 

Figure 1(a) is a schematic side-view and 1(b) is a top- view photograph of a circular kink 
formed on the sand surface. To measure the diameter of the kinks, a divider is placed near 
a kink and matched with the diameter. Several measurements are carried out for each kink 
and errors of measurements are about 1 mm. We also measure the surface shapes in some 
cases using the method of Ref. 6. 

First, we tilt the bucket by angle a and then turn on the motor to various speeds. 



In this type of experiment the initial granular surfaces are inclined flat surfaces. After a 
few minutes of rotation, we measure the radius of kinks r^. The radius depends on cu as 
rk oc (—0.2 < x < 0.1), thus a dimensionless radius Rk = VkUi^/g, which is roughly 

constant, can be introduced (Fig. 2). Next, the tilting angle is varied at fixed lo. When 
a < 30°, there forms a paraboloid- like granular surface whose shape depends on the initial 
condition but does not change with time. As a becomes larger than about 35°, a circular 
kink forms which is independent of the initial granular surface. Inside the kink, most of 
the grains on the surface avalanche during the rotation and we can see a complex motion 
like a whirlpool. The surface shape of the inner region is asymmetric about the rotation 
axis. On the other hand, those grains at and outside the kink are stationary with respect to 
the bucket. The surface shape of the solid-hke region is in general asymmetric and depends 
on the initial surface. However, when we first rotate the bucket vertically with an initially 
conical surface and then tilt slowly, the surface shape remains nearly symmetric. The kink 
is a boundary between the two dynamically distinguishable regions. As a increases, the Rk 
tends to increase. The measured values of Rk at several tilting angles are shown in Fig. 3. 
For a > 70°, a new type of instability appears; some sand grains are separated from the 
surface and fall freely during the rotation. 

We also study hysteresis by sequentially increasing, decreasing or randomly changing 
the angular velocity. For fixed a, the radius of kinks is determined only by cu and does not 
depend on the past history of changing cu. But the shape of the surface shows hysteresis. 
When uj increases a new kink appears inside the previous one. The previous kink formed at 
slower cu can be frozen in the solid-like outer region. Thus, there can be many concentric 
kinks. On the other hand, when cu decreases a new kink appears outside the previous one 
and the previous kinks are always washed out. 

Various vessels are tested as spinning buckets. If the width of a container is larger than 
the diameter of a kink, the radius of the kink does not depend on the size or shape of 
the container. We find that, to form the kinks, it is not necessary for the rotational axis 
to coincide with the axis of symmetry of a container. When we rotate a rectangular box 
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around an axis which is fixed on one of the side-walls, there appears a semicircular kink. 

Similar kinks are also found with the silica-gel (Matrex Silica, Amicon Corp., Danvers, 
MA 01923, U.S.A., of apparent density p = 0.35 g/cm^, particle diameter around 0.2 mm, 
and the angle of internal friction ^/ = 30 ± 1°), but the RkS are smaller by about 20 % than 
with the sand. We try sugar powder (sucrose, grain size ~ 0.2 mm, after 100°C for 2 hours 
to remove moisture), and observe kinks too. 

We now discuss the stability of the granular surface. There are four forces acting on 
a grain at the surface in the rotating frame: gravity mg, centrifugal force mfuj'^, normal 
force N and frictional (shearing) force f. Here, m is the mass of a grain, g the gravitational 
acceleration and r is the radial displacement of the grain from the rotational axis. We make 
the following assumptions. First, there is no bulk motion in the pile and the grain at the 
surface can only slide or roll on the surface. Second, there is no inertial effect. The grain 
stops as soon as it satisfies a force balance equation. Third, we also assume the Coulomb 
yield condition which states that the grain does not move when / < fiN, where n 



(/i = tan^j) is the coefficient of friction. The force balance equation for the granular system 
is 

N + mg + mfuj'^ + f = 0, f < fiN. (1) 
We define a dimensionless vector 

mg + mruj'^ _ ^ 

Uh = , (2) 

mg 

— * 

which is normal to the stable surface when there is no friction (/ = 0). The inclination 
of the rotational axis breaks the cylindrical symmetry of the vector flh and makes the angle 
P between fih and z time-dependent, where /3 is also the angle between the stable surface 
without friction and the bottom plane of the spinning bucket (Fig. 4). 

The kink is an abrupt change in the slope of the granular surface along the radial direc- 
tion. Therefore, we concentrate on the radial motion of grains. The radial stability condition 
can be expressed in terms of the angle f3 



taxiiPit) - 9f) < taxiO < tan(/3(i) + ^j), (3) 

where tan 9 is the local slope of the granular surface in the radial direction with respect to 
the bottom of the bucket (Fig. 4). 

The stable surface satisfying Eq. (3) for all t, can exist only if the following condition is 
satisfied; 

Pmax-Of < /^min + ^Z' (4) 

where /9max (/^min) maximum (minimum) value of (3{t) during rotation. If we use a 

cylindrical coordinate system (p, cf), z) aligned to the rotation axis, the and P{t) can be 
expressed as 

Hh = p(sin a sin — i?) + sin a cos + i cos a (5) 

and 

^/ X cosa . , 

V -R + 1 — 2Rsma sm cut 

where R = ruj^/g is a dimensionless radius and 6 = uot (the phase is chosen that = 7r/2 

corresponds to the highest position during rotation). The /?max (/^min) ('^) becomes 

the angle (5{t) at = 37r/2 (0 = 7r/2). 

Finally, we reach the following radial stability condition, 

R > \ I (7) 

_ Rc = ^ sin29/ («>^/)- 

If a < the steady state surface is stable for all i?. If a > the surface can be stable 
only in the region where R is larger than the critical radius Re- In the region with smaller R, 
avalanches of grains occur. We plot the Rc{a) in Fig. 3 for the sand sample with Of — 34° 
to compare with the radius of kink. The critical radii calculated by Eq. (7) reflect the 
qualitative features of the experimental kink radius; they increase with the inclination angle 
and do not depend on the angular velocity. 



The angular stability condition is also examined. The corresponding critical radius Rf for 
the azimuthal direction is always larger than Rc, and their difference increases monotonically 
as a increases. In our experiment (35° < a < 70°), however, the difference is small and the 
angular instability is expected to have little effect on the formation of kinks. 

In order to gain insights on the formation of the kinks, we study the system using a simple 



cellular automata model similar to that of Bak et al. |12| , p!3| . In addition to the assumptions 
made earlier, we further assume that the relative motion of grains in the azimuthal direction 
is not significant. Experimentally, grains just inside the kink show nearly circular trajectories 
without relative movement in the azimuthal direction. Since the kinks are formed by the 
motion of grains nearby the kinks, we expect that the above assumption does not alter their 
formations. The two-dimensional granular surfaces can now be described by one- dimensional 
curves — the surface profile h{R) at a given 0. 

The spatial coordinate R is made discrete — it is replaced hj Ri = i AR {i = —n, —n + 
1, ■ ■ ■ , ra), which runs from one end of the container to the other end. We measure local 
slopes Si, defined as {h{i) — h{i + 1))/AR, and check the stability of all the slopes following 
the criterion of Eq. (3). We then update the heights of the pile to make the local slopes at 
least marginally stable starting from the uphill and proceeding towards the downhill. 

If the local slope Sj_i becomes unstable due to an update at the i-th site, we allow 
"backward propagation" , where a perturbation at a site can influence a site uphill from the 
perturbation. To be more specific, we decrease the height h{i — 1) and transfer the excess 
amount to the {i + l)th site. If the change makes the {i — 2)th site unstable, its height 
is decreased in a similar way. We proceed until all local slopes become stable. Then the 
container is rotated by a small angle 6(f) and we again update the height. For computational 
simplicity we assume that the relaxation time of the pile is much shorter than the period of 
the rotation. At the boundary, we apply a mass conservation condition; grains cannot enter 
nor leave the container. 

In Fig. 5, we show time evolutions of the heights h{i,t) with different initial conditions 
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plotted during a half rotation. Here, AR = 0.01, a = 45°, Of = 34°, and 50 = f x lO^^ 
rad. One can see the region with large \R\ is solid-like; the profile does not change with 
time. On the other hand, the inner region is fluid-like; the heights keep changing. There is 
a discontinuity in the slope at the border of the two regions {\R\ ~ 1.2). The resulting kink 
looks similar to what is observed in the experiments. Although the formation of the kinks 
does not depend on the initial surface, the proflle of the solid-like region does. The radii 
Rk{oi) resulting from the simulation with Of — 34° are shown in Fig. 3 for different a in the 
range (0° — 90°). One notices that not only the quahtative features such as the threshold 
angle {a = Of) and the overall shape, but also their numerical values are in good agreement 
with the measured ones. However, this quantitative agreement may be accidental, since we 
expect the present simple model to reproduce qualitative, but not quantitative behaviors. 
To study the hysteresis behaviors, we start with a kink formed on the surface at a = 40°, and 
suddenly increase a to 65°. There indeed appears a new kink with smaller radius, leaving 
the old one stable, similar to what is observed in the experiments. 

Finally, the formation of the kinks can be viewed as follows. Consider the border at 
Rc from the stability analysis. There will be avalanches just inside the border. Due to the 
backward propagation, which we expect to exist in real sandpiles, the region just outside the 
border should be involved in the avalanche. This explains why the measured Rk is larger 
than Rc as shown in Fig. 3. The surface in the stable region will be marginally stable, so we 
expect its slope is a smooth function of R. The surface in the inner unstable region, however, 
is produced by a process different from that for the surface outside the border, so there is 
no reason to expect that the two slopes join smoothly. It is thus natural to expect that a 
sudden change in the slope is observed at the border between the sohd-hke and fluid-like 
regions. 

In summary, we have discovered circular kinks on the surface of granular material in a 
spinning bucket when its axis of rotation is tilted beyond the angle of internal friction. The 
radius of the kinks depends on the angular velocity, the tilting angle and the angle of internal 
friction. With a flxed tilting angle, the dimensionless radius of kinks Rk — TkUp' j g remains 
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roughly constant. We find that the surface is divided into two regions: a fluid-hke inner 
and a sohd-hke outer regions. We determine the critical radius Rc from the radial stability 
condition and the prediction reflects the basic features of the experiments. Using a simple 
cellular automata model, with the same stability condition and by allowing the propagation 
effect of avalanche, we obtain the granular surface in good accord with the experiments. 
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Fig. 1 




FIG. 1. (a) Schematic side-view of a tilted spinning bucket: the rotational axis z is tilted from 
the vertical direction by an angle a. A kink of radius rjt forms on the granular surface separating 
the solid-like (profile stable) outer region from the fluid-like (profile unstable) inner region, g 
denotes the gravity, (b) Top-view photograph of a kink formed on the sand surface (stressed by a 
dotted circle, of radius 2.0 cm, for a clearer view) rotated with uj = 230 rpm (direction is indicated 
by an arrow) and a = 45°. 
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Fig. 2 
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FIG. 2. The dimensionless radius = ^-^^ of the kinks formed on the sand pile. Examples 
for the two tilting angles, a = 60° (o) and a = 45° (•), indicate the relation Vk oc u!~^. Errors of 
the measurements are Ar^ 2± 1 mm. 
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Fig. 3 
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FIG. 3. The dimensionless radius Rk for sand at various a. Representative data for two values 

of u) are shown together with the critical radius i?c( ) calculated for stability and the Rk from 

simulation(— o — ) with the angle of internal friction (angle of repose) 9f = 34°. 
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Fig. 4 




FIG. 4. The vector fih is normal to the local surface (dashed line) when there is no friction. The 
angle f3 between the vector Hh and the rotation axis z is changing during rotation. The granular 
surface (solid line) in the shaded region is stable for the angle 6 between the granular surface and 
bottom of the bucket in the range (i — Of < 9 < (3 + 9 where ^/ is the angle of repose. This figure 
is an example that a grain is at the top position during rotation. 
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Fig. 5 




(b) 5 I 1 1 1 1 1 1 1 

-1.6 -1.2 -0.8 -0.4 0.0 0.4 0.8 1.2 1.6 

R 

Yoon et al. 

FIG. 5. Time evolutions of the granular surface calculated from the simulation of a simple 
cellular automata model with Of = 34° and a = 45° for two different initial conditions; (a) first 
tilt and then rotate, (b) first rotate vertically and then tilt slowly. Different initial preparations 
give different surface profiles, which however result in the same radii of kinks. The open and filled 
circles indicate the surface profiles at the initial and after a half (7r-radian) rotation positions, 
respectively, and two arrows represent the direction of changes. In (b), the filled squares are the 
experimental data, which are scaled by ^ , and the dashed line is the parabolic best fit to the data 
in the solid-like region (1.2 < \R\ < 4). The experimental data in the fluid-like region {\R\ < 1.2) 
are only rough estimates (Ah ~ 0.2) due to the steady flow of grains. The data at R > correspond 
to the upper half of the granular surface. 
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